FEB 



IN THE UNITED STATES PATENT AND TRADEMARK OFFICE 
re Patent Application of 
Junjie ZHU, et al. Date: January 25, 2002 

Serial No.: 09/883,683 Group Art Unit: 2631 

Filed: June 18, 2001 Examiner:--- 
For: METHOD AND APPARATUS FOR REDUCING SIGNAL DEGRADATION 



OFGS File No. : P/2778-13 ' JZ T 



Assistant Commissioner for Patents 
Washington, D.C. 20231 



SUBMISSION OF PRIORITY DOCUMENT 



Sir: 



009ZJd)U83^6ojOUM9ai 
2002 tZ a3d 
a3AI303d 



In accordance with 35 U.S.C. §119, Applicant confirms the prior request for 
priority under the International Convention and submits herewith the following document in 
support of the claim: 

Certified Singapore . Patent Application: 
(1) No. 200003402-5 filed June 16, 2000. 



I hereby certify that this correspondence is being 
deposited with the United States Postal Service with 
sufficient postage as first class mail in an envelope 
addressed to: Assistant Commissioner for Patents, 
Washington, D.C. 2023 1, on January 25, 2002: 

Robert C. Faber 
Name of applicant, assignee or 




January 25, 2002 
Date of Signature 



Respectfully submitted, 





Robert C. Faber 
Registration No.: 24,322 

OSTROLENK, FABER, GERB & SOFFEN, LLP 
1 180 Avenue of the Americas 
New York, New York 10036-8403 
Telephone: (212)382-0700 



RCF:sds 
Enclosure 



00524041.1 



REGISTRY OF PATENTS 
SINGAPORE 

This IS to certify that the annexed is a true copy of the foUowing 
Singapore patent application as filed in this Registry. 




Date of Filing 



16 JUNE 2000 



Application Niunber : 200003402-5 



Applicant(s) 



OKI TECHNO CENTRE (SINGAPORE) PTE 
LTD 



'tie of Invention 



: METHODS AND APPARATUS FOR 
REDUCING SIGNAL DEGRADATION 



RECBVED 

FEB 2 2 200^ 
Technologv Center 2600 



Sharmaine Wu Shee Mei 
Assistant Registrar 
for REGISTRAR OF PATENTS 
SINGAPORE 



200003A02-5 

1 6 JUN 2000 

Rule 19 



THE GRANT OF A PATENT IS REQUESTED BY THE UNDERSIGNED ON THE BASIS OF 

THE PRESENT APPLICATION. 



I, Title of Invention 


METHODS AND APPARATUS FOR REDUCING SIGNAL DEGRADATION 


II. Applicant (s) 

(See note 2) 


(a) Name 


Oki Techno Centre (Singapore) Pte Ltd 




Body Description / Residency 


A company organized and existing under the laws of 
Singapore 




Street Name & Number 


20 Science Park Road, #02-06/10, 
Teletech Park, Singapore Science Park II 




City 






State 






Country 


Singapore 1 1 7674 




(b) Name 






Body Description / Residency 






Street Name & Number 






City 






State 






Country 






(c) Name 






Body Description / Residency 






Street Name & Number 






City 






State 






Country 





PATENTS FORM 1 SINGAPORE 

THE PATENTS ACT 
(CHAPTER 221) 
PATENTS RULES 

The Registrar of Patents 
Registry of Patents 

REQUEST FOR THE GRANT OF A PATl 



2 00 003A02-5 



III. Declaration 


Country / Country Designated 




File No. 1 


of Priority 


Filing Date 








(see note 3) 


Country / Country Designated 




File No. 




Filing Date 










Country / Country Designated 




File No. 




Filing Date 










Country / Country Designated 




File No. 




Filing Date 









. Inventors 

{See note 4) 

(a) The applicant(s) is/are the 
sole/joint inventor{s). 

(b) A statement on Patents 
Fonn 8 is/will be furnished. 




□ YES 
H YES 



IE 
□ 



NO 
NO 



V. Name of Agent (if any) 
{See note 5) 


LLOYD WISE 


VI. Address for Service 

{See note 6) 


Block/Hse No 




Level No 




Unit No/PO Box 


P O BOX 636 


Postal 
Code 


910816 


Street Name 




Building Name 


TANJONG PAGAR POST OFFICE 


I. Claiming an earlier 
filing date under 
section 20 (3), 26(6) or 
47(4) 

{see note 7) 


Application No 




Filing Date 






[Please tick in the relevant space provided]: 

( ) Proceeding under rule 27( 1 )(a). 

Date on which the earlier application was amended = 

or 

( ) Proceeding under nile 27( 1 )(b). 



2 00 0(734 0 2-5 




VIII. Invention has been 
displayed at an 
International Exhibition 

(See note 8) 


□ YES E NO 


IX, Section 114 
requirements 

{See note 9) 


The invention relates to and/or used a micro-organism deposited for the purposes of 
disclosure in accordance with Section 114 with a depository authority under the 
Budapest Treaty 

□ YES HI NO 


X. Check List 

(To be filled in by 
applicant or agent) 


A. The application contains the following number of sheet(s) :- 


1 . Request. 


4 


Sheets 


2. Description. 


39 


Sheets 


3. Claim (s). 


9 


Sheets 


4. Drawing (s). 


6 


Sheets 


5. Abstract. 


1 


Sheets 


B. The application as filed is accompanied by :- 


1 . Priority document. 






2. Translation of priority document. 






3. Statement of Inventorship & right to grant 


X 




4. International Exhibition Certific^e. 






XI. Signature (s) 

{See note 10) 


Applicant (a) 




Date 


16 June 2000 




Applicant (b) 




Date 






Applicant (c) 






Date 







NOTES : 



2 00 005402-5 



I. This form when completed, should be brought or sent to the Registry of Patents together with the prescribed fee and 3 
copies of the description of the invention, and of any drawings. 

2 Enter the name and address of each applicant in the spaces provided at paragraph II. Names of individuals should be 

•licated in full and the surname or family name should also underlined. The names of all partners in a firm must be given in full. 
e place of residence of each individual should also be furnished in the space provided. Bodies corporate should be designated 
by their corporate name and country of incorporation and, where appropriate, the state of incorporation within that country should 
be entered where provided. Where more than three applicants are to be named, the names and address of the fourth and any further 
applicants should be given on a separate sheet attached to this Form together with the signature of each of these further applicants. 

3. The declaration of priority at paragraph III should state the date of the previous filing, and the country in which it was 
made and indicate the file number, if available. When the application relied upon in an Intemational Application or a regional 
patent application e.g. European patent application, one of the countries designated in that application [being one falling under the 
Patents (Convention Countries) Order] should be identified and the name of that country should be entered in the space provided. 

4. Where the applicant or applicants is/are the sole inventor or the joint inventors, paragraph IV should be completed by 
marking the *'YES' Box in the declaration (a) and the 'NO' Box in the alternative statement (b). Where this is not the case, the 'NO' 
Box in declaration (a) should be marked and a statement will be required to be filed on Patents Form 8. 

5. If the applicant has appointed an agent to act on his behalf, the agent's name should be indicated in the spaces available at 
paragraphs V. 

6. An address for service in Singapore to which all documents may be sent must be stated at paragraph VI. It is 
:ommended that a telephone number be provided if an agent is not appointed. 

7. When an application is made by virtue of section 20(3), 26(6) or 47(4), the appropriate section should be identified at 
paragraph Vll and the number of the earlier application or any patent granted thereon identified. Applicants proceeding under 
section 26(6) should identify which provision in mle 27 they are proceeding under. If the applicants are proceeding under rule 
27( 1 )(a), they should also indicate the date on which the earlier application was amended. 

8. Where the applicant wishes an earlier disclosure of the invention by him at an Intemational Exhibition to be disregarded 
in accordance with section 14(4)(c), then the 'YES' box at paragraph VIII should be marked. Otherwise the 'NO' box should be 
marked. 

9. Where in disclosing the invention the application refers to one or more micro-organisms deposited with a depository 
authority under the Budapest Treaty, then the 'YES' box at paragraph XI should be marked. Otherwise the 'NO' box should be 
marked. 

10. Attention is drawn to rules 90 and 105 of the Patent Rules. Where there are more than three applicants, see also Note 2 
above. 

I I . Applicants resident in Singapore are reminded that if the Registry of Patents considers that an application contains 
information the publication of which might be prejudicial to the defence of Singapore or the safety of the public, it may prohibit or 

itrict its publication or communication. Any person resident in Singapore and wishing to apply for paitent protection in other 
countries must first obtain permission from the Singapore Registry of Patents unless they have already applied for a patent for the 
same invention in Singapore. In the latter case, no application should be made overseas until at least two months after the 
application has been filed in Singapore. 



For Official Use 



Application Filing Date 
Request received on 
Fee received on 
Amount 

* Cash / Cheque / Money Order No : 



/ / 
/ / 
/ / 



* Delete whichever is inapplicable. 



Patents Form 1 —4 



Our Ref: SP3899 



SP389904.doc 



2 l5 

2 00 0 0 3 4 0^-i 

1 6 JUN 2600 



Methods and Apparatus for Reducing Signal Degradation 



Field of the invention 

The present invention relates to methods and devices for reducing noise in a 
5 signal received and re-transmitted by a relay station. The invention is 
particularly suitable if the signal received by the relay station includes 
multipath interference. The invention further relates to a relay station, such as 
a single frequency relay station. 

10 Background 



The ISDB-T (Integrated Services Digital Broadcasting - Terrestrial) BST- 
OFDM (band segmented transmission - orthogonal frequency division 
multiplexing) system is an example of what is often called a "single frequency 
15 network" (SFN). 



In fact, the signal is generally made up of multiple carrier signals. Each 
carrier has a finite, though narrow, frequency range, and when we talk about 
the frequency of a single carrier signal, we mean the center frequency of the 
20 narrow range. For multi-carrier transmission, such as OFDM in an ISDB-T 
system, each carrier has a different respective center frequency, and the 
frequencies are orthogonal to each other The multi-carrier signal is the sum 
of the signals of the multiple carriers and is a wideband signal which has a 
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center frequency too. This center frequency is here referred to as the 
frequency of the signal. 



In conventional broadcasting, a signal received by a relay station at one 
5 frequency is re-transmitted at another frequency in order to avoid 
interference. However, by contrast, in digital broadcasting, because of the 
scarceness of frequency resources, it is usually required that sanne signal is 
re-transmitted at the same frequency (i.e. over the same range of frequency) 
by all transmitters. This is what is meant by the term "single frequency 
10 network". 

Figure 1 shows a block diagram of an example of the configuration of a 
known relay station (also referred to here as a repeater). To relay the 
terrestrial digital broadcast with SFN, waves 1 received from a master station 
15 by an antenna 3 are amplified by a relay broadcasting system 4 and re- 
transmitted by a second antenna 7 at the same frequency. The waves 1 
include a master signal and a level of distortion. 

The field strength of the transmitted wave is very intense compared with that 
20 of the received wave. Although the receiving antenna 3 is directed to the 
master station, the transmitting antenna 7 is arranged with a directivity which 
matches the shape of its service area, and the radio wave is normally 
radiated at a wider angle. Therefore, coupling waves 5 are received by the 
receiving antenna 3 at a level which increases with the directivity index of the 



transmitting antenna 7 to the direction of the receiving antenna 3, or 
conversely with the directivity index of the receiving antenna 3 to the direction 
of the transmitting antenna 7. Furthermore, the coupling waves 5 include 
components of the signal transmitted by the transmitting antenna 7 which 
reaches the receiving antenna 3 after being reflected by structures, trees and 
mountains near the relay station. Consequently, the coupling waves 5 are a 
synthetic signal having a varying delay time, intensity and phase. 

A repeater including a cancellation algorithm has been proposed. As shown 
in Figure 1, the amplification system 4 of the relay station includes not only 
an amplifier 13 but also a coupling waves cancellation unit 9 and a BPF 
(bandpass filter) 11. Note that the apparatus in Figure 1 is simplified, and in 
practice often includes additional components such as D/A converter(s). 

The coupling waves cancellation unit 9 is convenfionally designed based on 
the assumpfion that the distortion of the incoming radio waves from the 
master stafion is small compared to the master wave itself. However, this 
assumption is not appropriate for a practical repeater transmission 
environment which generally has a mulfipath transmission channel (e.g. there 
are multiple radio paths between the master station and the relay station, 
each taking different times). In a mulfipath channel, the distortion or 
attenuafion of the master wave can be serious at a particular frequency (i.e. 
over a particular frequency range), especially when the amplitudes of the 
paths arriving at different fime are the same. This is called frequency 




selective fading. In the presence of deep frequency selective fading caused 
by such multipath interference, the conventional coupling cancellation unit 
cannot work well and diverges, because the conventional coupling 
cancellation unit algorithm does not take the effect of multipath interference 
5 into account. Furthermore, it is very difficult for a relay station to distinguish 
between multipath interference and coupling signals in the received signals. 
Moreover, for non-coherent modulation, the conventional PSK (phase shift 
keying) estimation method cannot work well if there is a multipath channel, 

10 Both coupling interference and multipath interference vary with time. Thus, a 
further advantageous feature of a cancellation unit is that the hardware 
processing speed should be fast enough to update the adaptive filter of the 
coupling cancellation unit. The time taken depends on the amount of 
computation required to perform the updating and the hardware available. 

15 

Summary of the present invention 

The present invention aims to provide a new and useful method of operating 
a relay station, and a cancellation circuit for a relay station. 

20 The present invention proposes in general terms that the transmitted signal is 
interrupted (e.g. by interrupting the connection between the cancellation 
circuit and the amplifier) so that the multipath interference can be estimated 
in the absence of coupling. Subsequently, the transmission is resumed, and 
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the digital filter is adjusted to take into account both the multipath interference 
and the coupling. 



Specifically, in a first aspect the invention proposes a method of nriodifying 
5 the values of a plurality of digital filter coefficients for use by a digital filter 
which is a component of a relay station, the relay station successively 
receiving a first signal, modifying the first signal using the digital filter to form 
a second signal, and transmitting the second signal with amplification, said 
first signal including a master signal subject to an interference and a coupling 
10 component derived from the second signal, the method comprising: 

(i) interrupting said transmission of the second signal; 

(ii) deriving an estimate of said interference and storing the derived 
interference estimate; 

(iii) resuming said transmission of the second signal; and 

15 (iv) modifying said digital filter coefficients using said stored 

interference estimate and said second signal, to reduce the difference 
between said second signal and said master signal. 

Step (iv) may be repeated iteratively so that the difference between said 
20 second signal and said master signal is successively reduced. To monitor the 
degree to which the relay circuit is successful in reducing the (e.g. multipath) 
interference, a first parameter of the relay system characterizing the 




difference between said second signal and said master signal is repeatedly 
derived. Step (iv) may be repeated until the value of said first parameter 
traverses a predefined convergence value. Thereafter adjustment of the filter 
coefficients may be stopped, or optionally further adjustments may be carried 
5 out, e.g. to compensate for variations in the coupling component, neglecting 
the stored interference estimate. This may be adequate while the multipath 
interference remains approximately constant, but as the multipath 
interference varies it will become less adequate. Accordingly, a second 
parameter of the relay station, also characterizing the difference between the 

10 second signal and the master signal, may be repeatedly derived, and when 
the second parameter traverses a predefined divergence value, said steps (i) 
to (iv) may be repeated. Note that the first and second system parameters 
may be the same, but in this case the divergence and convergence values 
are chosen to be different, so that there is a range of values for the system 

15 parameter in which steps (i) to (iv) are not being carried out. 

Preferably, the relay station modifies the first signal by passing it through a 
coupling cancellation unit. This has an input receiving the first signal and an 
output outputting said second signal. The cancellation unit includes the digital 
20 filter and a subtracter. The input of the coupling cancellation unit is an input 
to the subtracter and an output of the subtracter is passed to the output of the 
cancellation unit and to an input of the digital filter. The subtracter subtracts 
the output of the digital filter from the received signal. 
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The relay station further includes an amplifier to amplify the second signal 
before it is transmitted. A switch may be provided between the cancellation 
unit and the amplifier, as the means for interrupting the transmission of the 
5 second signal. 

The step of modifying said digital filter coefficients may be performed, as 
explained in detail below, by operations in frequency space which produce 
results which are then converted into time space by a discrete inverse 
10 Fourier transform, and then used to modify the coefficients. 

To adequately remove muitipath interference, the number of digital 
coefficients (the tap length, M) itself need be no higher than approximately 
the product of a first value (e.g. ^5yls) representing the maximum delay of the 
15 muitipath signal and the FFT sampling frequency (fs). 

Thus, whereas conventionally the tap weights are found by an interpolation 
over all the carriers of the signal, the present inventors have recognized that 
this is inefficient. A much smaller number of virtual carrier can be used. 
20 Specifically, the inverse Fourier transform according to the present invention 
preferably generates a number of time domain points which equal to or a little 
larger than (e.g. up to twice or three times) the tap length of the digital filter. 
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To generate this number of output time domain points to sufficient accuracy, 
the IFFT need only operate on a number of frequency domain points ("virtual 
carriers") which is of the same order, such as the number of output time 
domain points multiplied by the ratio of the number of carriers to the FFT 
5 size. 

This produces a significant improvement in efficiency in relation to the 
conventional method, and in effect constitutes an independent second aspect 
of the invention, which is applicable even in cases in which transmission of 
10 the second signal is not interrupted as explained above. 

According to this second aspect, the invention proposes a method of 
modifying the values of a plurality of digital filter coefficients for use by a 
digital filter which is a component of a relay station, the relay station 

15 successively receiving a first signal, modifying the first signal using the digital 
filter to form a second signal, and transmitting the second signal with 
amplification, said first signal including a master signal subject to an 
interference having a maximum delay time, each of the digital coefficients of 
the digital filter modifying the first signal at pairwise time delay intervals Tp, 

20 the method including: 

deriving a discrete Fourier transform of the second signal; 
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employing said discrete Fourier transform of the transfer function to 
derive correction values for said digital filter coefficients expressed in the 
frequency domain; and 

transforming said correction values into the time domain by an inverse 
5 Fourier transform operation, the digital filter coefficients being modified based 
on the time domain correction values, 

wherein: 

the product of Tp and the number of digital coefficients is substantially 
equal to the maximum delay time of the interference component of the first 
10 signal, and 

the inverse Fourier transform generates an number of output time 
domain values equal to or slightly larger than the number of digital 
coefficients. 



15 In either aspect of the invention, preferably, the master signal is in a format 
(differential modulation) in which symbols are coded by differences between 
consecutive segments of the master signal. The adjustment of the digital filter 
coefficients may then be performed by a process which is based on 
estimating the differences between consecutive segments of the master 

20 signal, i.e. neglecting an estimation of what the master signal was in absolute 
terms. For example, the master signal may be in a PSK format in which the 
master signal encodes symbols by phase differences between consecutive 
segments, and the method may include estimating the differences between 
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the phases of consecutive segments, e.g. rather than trying to estimate the 
phase of each segment individually. This is particularly useful in determining 
the transfer function of the relay station. 

In fact, this technique constitutes an independent third aspect of the 
invention, which is a method of modifying the values of a plurality of digital 
filter coefficients for use by a digital filter which is a component of a relay 
station, the relay station successively receiving a first signal, modifying the 
first signal using the digital filter to form a second signal, and transmitting the 
second signal with amplification, said first signal including a master signal 
and a coupling component, the master signal being in a format in which 
symbols are coded by differences between consecutive segments of the 
master signal, and the coupling component being derived from the second 
signal, the method comprising: 

(i) deriving an estimate of differences between consecutive 
segments of the master signal; 

(ii) deriving an estimated error in said estimation of differences 
between consecutive segments of the master signal; and 

(iii) modifying said digital filter coefficients to minimize said 
estimated error. 

The invention also relates to a cancellation unit arranged as discussed 
above, as well as to a relay station incorporating such a cancellation unit. 
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Brief description of the figures 

Embodiments of the invention will now be explained for the purpose of 
example only, with reference to the following figures, in which: 

Figure 1 is a block diagram of a relay station system; 

Figure 2 is a block diagram of a relay station including a coupling 
cancellation unit according to an embodiment of the invention; 

Figure 3 is a block diagram of using digital signal processing to 
estimate FIR filter coefficients in the arrangement of Figure 2; 

Figure 4 illustrates a method for calculating the transfer function for 
PSK-OFDM; 

Figure 5 shows the layout of a pilot signal; 

Figure 6 shows simulated learning curves for a 64QAM-OFDM 
cancellation unit, using IFFT algorithms with 2048 points and 256 points 
respectively, for DUm=6dB, DU=2dB, |.i=0.5 and two values for the coupling 
delay; 

Figure 7 shows simulated learning curves for a 64QAM-OFDM 
cancellation unit, with and without channel pre-estimation, for DUm=9dB, 
DU=2dB, 1.1=0.5, and IFFT of 256 points; 
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Figure 8 shows simulated learning curves for a 64QAM-OFDM 
cancellation unit with and without channel pre-estimation for DUm=5dB, 
DU=2dB, |a=0.5, and IFFT of 256 points; 

Figure 9 shows simulated learning curves of a 64QAM-OFDM 
5 cancellation unit in a multipath channel with deep fading and channel pre- 
estimation; 

Figure 10 shows the definition of the SN ratio at a receiver; 

Figure 11 shows simulated learning curves of a tc/4-DQPSK-OFDM 
cancellation unit in a multipath channel with various algorithms for DUm=7dB, 
!0 DU=2dB, ).i=1 , and IFFT of 256 points; 

Figure 12 shows simulated learning curves of a tiM-DQPSK-OFDM 
cancellation unit in a multipath channel with various algorithms for DU=2dB, 
).i=1 , and IFFT of 256 points; and 

Figure 13 shows simulated learning curves of a k/4-DQPSK-OFDM 
15 cancellation unit in a multipath channel with various algorithms for DU=2dB, 
1-1=0.5, and IFFT of 256 points. 



Detailed description of embodiments 

Like the conventional system discussed above, the present embodiment has 
20 the general structure illustrated in Figure 1, but the coupling waves 
cancellation unit 9 is as illustrated in Figure 2. The optional BPF 11 of Fig. 1 
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is omitted from Fig. 2 for simplicity. Likewise, Figure 2 omits any A/D 
converter(s) in the system. 



Figure 2 shows the detail of the coupling wave cancellation unit and the 
5 principle of canceling coupling waves between the transmitting and receiving 
antennas of the broadcasting wave relay station. The master signal is x(t) (in 
the frequency domain X(co)), which is subject to multipath interference h(t) 
(H(co), in the frequency domain, is the multipath channel transfer function), in 
the channel 15 between the master station and the relay station, to produce a 
10 signal r(t) (R(co) in the frequency domain). 

The signal R(co) from the master station is given by 

R{co)=^X{o))h{co) (1) 

A coupling wave is added to r(t) (a process represented here by the notional 
adder 17), to produce the signal received by the antenna 3, here called i(t) 
15 (l(co) in the frequency domain). 

A subtracter 19 is used to subtract the coupling signals and multipath 
interference from the incoming signal i(t) which is one input of the subtracter 
19. An adaptive FIR filter 21 is fed the output of the subtracter 19 and 
20 connected to the subtracting terminal of the subtracter 19. And a digital signal 
processing block 23, which estimates the characteristics of coupling path and 
multipath, is fed the output of the subtracter 19 and provides the FIR filter 21 
with updated FIR filter coefficients. Moreover, a switch 25 is inserted between 
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the subtracter 19 and the amplifier 13 to control the relay of the broadcasting 
signal. When the switch 25 is switched on, connecting the subtracter 19 to 
the amplifier 13, the signal after cancellation through the subtracter 19 is 
transmitted through the antenna 7. When the switch 25 is switched off, 
switching off the connection between the subtracter 19 and the amplifier 13, 
the output signal of the subtracter 19 is not transmitted and hence there is no 
coupling wave in the incoming signals. 



Supposing the characteristics of coupling wave frequency are C(o)), the gains 
10 and frequency characteristics of amplifier at relay station 13 are G(o)). 
Because the system is studied in baseband, G(co) may be treated as not 
including an effect from the antenna or any up/down converter. The signal 
frequency characteristics at observation point P are written as S(co), the 
incoming signal at the relay station receiving antenna 3 is /(co), expressed by 
15 the following equation: 

j{co) = S{co) G{a}yc{co)^ R{co) (2) 

Now supposing that l/l/(a)) is the frequency characteristics of the cancellation 
unit, 

S{co)^l{co)-S{coyw{co) (3) 
20 Substitute equation (2) for (3). and rearrange it as follows: 

^fe) (4) 
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By further substituting equation (1) for (4), we express the equation (4) as 
follows: 



5(.) = __^M_£M^ (5) 
^ ' \-{G{(jo\C{(o)-W{co)] ^ ' 

Consequently, transmission function F(co) of whole the system is expressed 
5 with the following equation: 

F{co) = 4^ = , , , (6) 

^ ^ X{co) \-{G{o)yc{(D)-W{co)] 



Using this formalism, we can analyse the behaviour of the conventional 
method of canceling coupling waves discussed above. In the conventional 
10 method it is assumed that the distortion oftfie incoming radio waves from tfie 
master station is small enough. That means that H(co) = 1 . Then the system 
transfer function can be written as 

In this case, to cancel the coupling waves, G{co) c{co)=lv{co) is required. In 
15 other words, the optimal FIR filter transfer function is IV^^^X^^)^ G{^)'C{co) . 
Now setting the tap-weight estimation error at the /+1-th iteration is 

Er{co)=W^JcoywM (8) 

where f^^+,(6;) = ff'^,(^y)+ A^K.^,(rf;) is the updating tap-weight vector at /-th 
iteration. Obviously, if the change in the tap-weight vector at /-th iteration 
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A^F;.^i(^y)= and Er{oi) is controlled to the minimunfi, coupling can be 

cancelled. It is expressed by the following equation, which can also be 
derived using the NLMS (Normalized Least Mean Square) algorithm (see 
Simon Haykin, "Adaptive Filter Theory", Third edition, 1996): 

By introducing a scaling factor |a, we may write 

PK,,,M^fK,M + //-APF,,,(^) (10) 

This method is limited to the characteristics of multipath channel. We note 
that in Equation (9), F(co) should not be very small, or the cancellation unit will 

10 become very unstable and hence diverge. From Equation (6), we know that 
F(a)) is a function of multipath channel, /-/(co). In multipath channel, when the 
attenuation difference between two main paths is not very large (that is, 
"DUm", defined as the power ratio of the first path to the second path, is 
small) there is a deep fading at some frequency and /-/(o)) at that frequency 

15 will be very small. Therefore, it is possible that the value of F(co) at some 
frequency is very small if /-/(co) is small. 



Let us now compare this with the algorithm proposed by the present 
invention. To do this we analyse Fig. 2 again in the context of the present 
20 invention. As we know, the target of the coupling cancellation unit 9 is to 
cancel the coupling waves and the effect of multipath channel and hence 
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recover the original signal from nnaster station. That is, the output of coupling 
cancellation unit 9, S(a)), will converge to be almost equal to master station 
output. X(co), i.e. 

F(a)) = 1 (11) 

To hold this condition, the optimal FIR filter transfer function is developed as 
following. 

Kr, (^) = G{cd)c{co)-^ H{co)-\ (12) 

Using Equation (6) and (12) in (8), we may reformulate the estimation error 
Er{co) as 

Er{o}) = G{co)^ C{o)) + H{o)) - 1 - AW.^, {co)~ W. {co) 

If Er(co) is controlled to the minimum, coupling can be cancelled. Then we 
hope that APF.^,(6>)is minimized, which is subject to the constraint 

Er{co) = 0 (14) 
Then Equation (13) is converted as follows: 



(15) 



This result is a little different from that of conventional method in which the 
multipath channel effect, /-/(co), is not considered. Now let's see what will 
happen if H((o) is ignored. 
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X{co) 



F{co) ~ S{co) 
H{co)^G{co)' C{co)-\ - W{co) 
H{co) 



(16) 



From the equation (16), we know that when H/((o) approaches WopM, 
^^^;4-i(^) approaches zero. In this case, the multipath channel, H(to), has little 
influence on A(F.^,(a>). Now it can be concluded that in multipath channel 
5 when H(o)) approaches 1. the coupling cancellation unit 9 can work perfectly; 
but when the channel has large frequency selective fading, the cancellation 
unit 9 may be very unstable. 



From Equation (15), we know that if we can estimate the multipath channel 
10 information, we can apply the channel function /-/(co) to cancellation unit 9 and 
improve the system performance. However, it is very difficult to implement 
the equation (15) directly. Since for repeater the received signal includes not 
only the coupling wave but also the multipath interference, it is almost 
impossible to distinguish the two factors and extract /-/(w). Therefore, the 
15 problem is how to obtain the multipath channel information. 



The embodiment of the present invention is operated by switching off the 
switch 25 to stop the connection of the coupling cancellation unit 9 and the 
amplifier 13 in Figure 2. Thus, there is no signal re-transmitted at the relay 
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station and hence no coupling waves in the inconning signals at the relay 
station. The estinnation of multipath channel can be processed. 



The new proposed method is called Channel Pre-estimation (CPE): we 
5 control the repeater not to transmit signal until the multipath channel is 
estimated. In the repeater environment, it is usually assumed that the 
transmission channel is stationary for a long time and there exists a line-of- 
sight (LOS) from the master station antenna to the relay station antenna. This 
implies that multipath DUm is larger than zero and the channel /-/(co) is 
10 constant for a long time within which cancellation unit can converge 
completely. At the beginning, the repeater doesn't retransmit the signal from 
the master station. Since at this moment there is no coupling, the repeater 
can estimate the multipath channel information, H(co). The required duration 
of time to stop transmission to estimated multipath channel information is 
15 only several milliseconds. After starting to repeat the incoming signal, the 
coupling cancellation unit can use /-/(co) to compute the tap-weight W{io) with 
Equation (15). 

The digital signal processing block to implement the CPE method is shown in 
20 detail schematically in Figure 3, wherein the incoming signals (at point P) are 
transformed into frequency domain using a FFT (step 29). For QAM 
modulation, for example, the signals over pilot carriers are extracted (step 
31) for estimation (by interpolation in step 33) of the channel transfer function 
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(CTF) F. There is a switch 35 connecting the interpolated CTF to the 
multipath channel estimation or the FIR filter coefficients updating estimation. 
The operation of the switch 35 in Figure 3 corresponds to that of the switch 
25 in Figure 2. When output of interpolation block 33 is connected to 
5 multipath channel estimation block 37, the multipath channel estimation is 
processed and multipath channel information /-/(co) is stored. When the 
output of interpolation block 33 is connected to the filter coefficients update 
estimation block 39 which is also fed with the estimated multipath channel 
information, the FIR filter coefficents update estimation is processed. Then 
10 the updates to the filter coefficients are transformed into time domain using 
IFFT (step 40) and a correction value is extracted (step 41), multiplied by a 
coefficient (step 42), to produce a correction value which is used for 
correcting the FIR filter values (step 43). 

15 Thus, the method of canceling coupling waves includes the steps of: 

a) switch off the connection 25 of the coupling cancellation unit 9 and the 
amplifier 13 to make sure that there is no signal relayed from the 
repeater; 

b) when there is no signal repeated, estimate the multipath channel H(co) 
20 and store the estimated channel information; 

c) Switch on the connection 25 of the coupling cancellation unit 9 and the 
amplifier 13 to allow the output signals of the coupling cancellation unit 9 
to be relayed; 
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d) apply the estimated multipath channel information to estimate the FIR 
filter coefficients; 

e) Update the FIR filter coefficients; 

Rather than a switch 25, the relay station may employ any other mechanism 
to disable the amplifier from transmitting an amplified signal to the antenna 7 
while the filter coefficients are being tested. 

The steps of estimation of FIR filter coefficients can be repeated at any time. 
We may arrange a device to measure some parameter, e.g. Signal-to-Noise 
(SNR) of the coupling cancellation unit output or mean square error of FIR 
filter coefficients. If the parameter rises to a divergence threshold, e.g. 
because the multipath channel has changed since the FIR coefficients were 
last derived, this indicates that the coupling cancellation unit 9 may diverge. 
In this case, the repeater stops retransmission, estimates the multipath 
channel again, and makes the cancellation unit 9 adapt to the channel and 
converge to the required extent. 

As we know, after convergence when l/l/(co) approaches l/l/op/(o)), AlV^^^^) 
approaches zero. At this moment, the multipath channel, /-/(co), has little 
influence on the estimation of AW^^X^). Then, the estimation of multipath 
information is not needed and a function for the channel information, H(co), 
can be assumed, e.g. it can be assumed to be one. 
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A convergence threshold of some system parameter (e.g. the same 
parameter which triggers the recalibration process by reaching the 
divergence threshold) is used to distinguish whether the coupling 
cancellation unit 9 has yet reached this state of convergence. When the 
5 value of the measured parameter of the coupling cancellation unit 9 reaches 
the required threshold, the operation discussed above (i.e. neglecting further 
estimates of H(co)) can be executed. Note that the value of the convergence 
threshold may be different from that of the divergence threshold. 



10 From Equation (15), we can see that if a judged value of the pilot signal or . 
received signal is used to determine X(o)) by measuring the frequency 
characteristics S((o) of the signal at the observation point P, the coupling and 
multipath interference can be cancelled. In ISDB-T systems, several 
modulations are used, nf4 Shift DQPSK, QPSK, 16QAM and 64QAM. For 

15 non-coherent and coherent modulation, the methods of estimating X{w) used 
in the present invention could be different. We will now discuss these 
methods in detail for the cases of differential modulation, and modulation with 
scattered pilot signals. 
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1. Differential Modulation 

For differential modulation such as tc/4 Shift DQPSK, the desired signal. X(to), 
can be.estinnated using conventional PSK estimation method as shown in 
Figure 4. 

5 The phase and amplitude for the /c-th carrier of /-th received symbol can be 
obtained from Equation (17) and (18). 

0Ai,k) = tan- ^gi^ (17) 

r(i, k) = 4si\i.k)-^SQ\i,k) ( 1 8) 

Here SI and SQ are defined as the quadrature components of S, i.e. the real 
10 and imaginary parts of the complex signal S respectively. The phase 
difference between the received signal and the estimated desired signal is 
given by 



e{i,k)^eXhk)-^RND 



M,-eXi,k) 



2n 



(19) 



where RND[X] rounds the elements of X to the nearest integers. The 
15 amplitude of desired signal as unit one. Then system transfer function 
defined by Equation (6) is written as 

F{i,k) = r(i,k)txp{}0{i,k)) (20) 

As we know, the positions of k/4-DQPSK signals in the complex coordinates 
alternate between constellation Podd = {-3n/4, StiM. -tiM, ttM} and Peven = {0, 
20 71, 7i/2, -7r/2}, In practical systems, the receiver doesn't know which 
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constellation the received signal belongs to. Therefore, an eight-point 
constellation should be used in PSK-estimation. i.e. Mp = 8 in Equation (19). 
This results in higher estimation error probability and deteriorates system 
performance. 

Here, we just consider the phase difference between the consecutive 
DQPSK-OFDM symbols, 

A0(i + ],k) = ii+i,k)- (/, k) (21) 

Because the original DQPSK data and its demodulation are related to the 
phase difference only, we do not have to estimate the phase of every 
received symbol like in PSK estimation method. Then we may estimate the 
desired differential phase. A<9(/ + 1,^), by using the same algorithm as PSK 
estimation for individual symbol phase, 



M{i+\,k) = — -RND 
M,. 



A/,.A(9(/ + 1,A-) 
In 



(22) 



In order to recover the original data, we consider the estimation error of the 
phase difference, 

he, (/■ + I, /r) = A(9(/ + 1. A:) - A<9(/ + 1, k) 

= (a^(/ + 1, At) + ^, (/, k)) - (A9{i + \,k) + 0, {i,k)) (23) 
= <?, (/ + 1, *) - (a^(,- + ],k) + 0, (/, k)) 

We define the amplitude of DQPSK signal is unit one and the phase of the 
estimated desired signal is 

0y ii+lk) = e{i,k) + A^(/ + 1, A:) (24) 
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then the estimated desired signal is written as 

X(i + lk) = exp(^,, (/ -f 1,/:))= exp(0{i,k) + A0(i + 1,^) (25) 

Substitute Equation (24) into (23), 

A0^ (/■ + I, /:) = (/ + I, A:) - 0, (/ + 1, *) (26) 

5 Therefore, the target of differential nnodulation coupling cancellation unit is to 
make A6^^(/ + l,/:) to the minimum. This method can be implemented by 
replacing 0,{i,k) and 0(1, k) with A0(i + \,k) and A<9,(/ + l,A:)in Figure 4. 

Note that the estimated desired signal may not be equal to the original signal 
which is transmitted by master station and distorted by multipath channel. 
10 The system transfer function is represented as 

Xii,k) 

= rii, k) exp(j(^,. (/, k) - e, (/, k))) 

= rii, k) exp[j(^., (/, k) - 0,, (/, k)\ exp[j(^., (/, k) - 0, (/, ^)) (27) 
= /•(/,/:)exp(jA^,(/,A:)) exp[j((9.,(/,A) 
= • exp[j(^,. (/, k) - e, (/, -fe))] 

Obviously, it is impossible to calculate F{i,k), since we do not know the 

original signal phase, 6^{i,k). However, we can still use F\i,k) = ^M^ by 

X(,i,k) 

replacing F{i,k) \\\ Equation (15). 
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This method is better than the conventional method. Because in a multipath 
channel, the transmitted signal can be seriously distorted, the error 
probability of the PSK estimation of desired signal would be large especially 
at Mp = 8. But assuming the multipath channel is quasi-static during two 
5 symbols, there is little influence on the phase difference between the 
consecutive symbols. Another advantage is that the phase difference 
between consecutive k/4-DQPSK signals is always from a four-point 
constellation P= {-37t/4, SttM, -kM, uIA). In this case, the minimum distance 
between the pairs of signal points is enlarged, i.e. Mp - 4. Therefore, the 
10 differential PSK estimation is more reliable than the PSK estimation. 

However, in a worse multipath channel, neither PSK nor DQPSK estimation 
can really realize the transmitted signal due to lack of the knowledge of the 
real initial phase of the transmitted signal. Therefore, the ability to equalize or 
15 cancel the multipath interference for the cancellation unit is limited by the 
depth of multipath fading. If the initial phase of the transmitted signal is 
required, a training sequence for estimation of the initial phase should be 
inserted before transmitted data. 

20 2. QAM Modulation with Scattered Pilot Symbols 

In ISDB-T systems, for coherent modulation segments, the SP (Scattered 
Pilot) signals (i.e. data symbols with an amplitude value known at the 
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receiver) are transmitted together with data carriers, so that the channel 
transfer function (CTF) is estinnated both in frequency and in time. Figure 5 
shows the scattering pattern of pilot carriers distributed in the direction of 
OFDM symbol and carrier. According to the OFDM transmission system, pilot 

5 and data are regarded as signals under the frequency domain. Transmission 
is made after making IFFT on the modulation side and FFT on the 
demodulation side. Since the transmission signal X(/, k) is known in the pilot 
signal section, the samples of the CTF, F{i\ k) can be easily obtained by 
dividing the received pilot signals by the known pilots. The channel 

10 estimation over all carriers is performed by interpolation of the CTF samples. 
The interpolation processing is made using two-dimensional filters after 
estimating the pilot signal as a reference. 



After obtaining t:^W-^^[co), IFFT and coefficients extraction is performed 

1 5 Am</) = Coefficients - Exfraction{lFFT[AJV.^, (co)]} (28) 

where Aw(i) is a (/W x 1) vector, and M is tap length of the FIR filter. Since 
there is a FFT sample delay between the input and output of FIR filter, Awji) 
coefficients are extracted from the second samples to the following M+1 
samples after IFFT. Finally, we get 

20 M</ + 1) = w(/) + //Au</) (29) 
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The IFFT used in any of nnethods according to the present invention 
discussed above may be a conventional IFFT algorithm, examples of which 
5 are well- known in this field. However, we will now explain why it is possible 
to improve on the efficiency of such a strategy. 

In the coupling wave canceller, an adaptive filter characterizing the coupling 
path and multipath channel is used. Two important parameters of the 

10 adaptive filter are tap length M and tap interval Tp. The tap length is the 
number of taps in the filter. The tap interval refers to the delay time between 
the consecutive taps of the filter. The two parameters determine how long a 
delay time of the multipath signal or coupling signal the filter can 
accommodate. In other words, the tap length should be longer than the 

15 maximum delay time. 

In a repeater environment, for example, the maximum delay time of the 
multipath channel is normally up to about 15|.is. To cope with such a 
multipath channel, M Tp should be larger than 15|_is. Because FFT (2048) 
20 sampling frequency fs in ISDB-T is 8.1 2693. ..MHz, the sampling interval Tc = 
1/ fs and the FIR filter tap interval should be Tc too. Thus, M > 15i_is/ Tc and 
tap length is set to 128. 
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The adaptive filter coefficients are generated through a series of connputation 
over frequency donnain and transforming the computation results into time 
domain using IFFT. As we know, the IFFT output sample point number is a 
number of the n-th power of 2. It could be very large when its input point 
number very large. When the output IFFT size is 2048 points, only the first 
128 IFFT output samples out of 2048 need to extracted for use as the filter 
coefficients. 

Now let us consider the conventional technique of making use of all of the 
IFFT output points for the filter coefficients. Then the tap length is equal to 
the IFFT size 2048. Using the example above, such an adaptive filter can 
adapt to a multipath delay time of 20487c = 252|.ls long. This is surplus to 
requirements, because the multipath signals received after 15|.is are 
insignificant, and can be ignored. In other words, we do not need to compute 
an IFFT with a size up to the conventional value 2048. It is enough to use a 
filter with 128 taps to obtain only previous 15|.ls channel information. 
Therefore, we may reduce the IFFT size to a value equal or a little larger than 
the tap length. 

Theoretically, it is possible to exactly reconstruct a band-limited function if the 
sample interval ts is smaller than 1/2fs where fs is the frequency limit of the 
function, but ideal interpolation must be used. In ISDB-T channel estimation, 
since the interpolation is done in the frequency domain, the sample interval 
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for the interpolation of the sampled CTF is the spacing between the pilot 
carriers, Af, while the frequency limit is the coherent bandwidth of the channel 
Be. In all practical channels the coherence bandwidth of the channel Be is 
larger than the pilot spacing Af thus the CTF may be considered as 
oversampled and can be restored over all carriers. Furthermore, the channel 
impulse response can be obtained after the IFFT transform of the 
interpolated CTF. 

In the coupling cancellation unit 9. the tap coefficients of the FIR filter are 
extracted from the IFFT transform of the estimated CTF. The filter tap-length 
is determined by the delay time of the coupling waves which is related to the 
coherent bandwidth of the transmission channel. Because normally tap 
length, e.g. 128, is much smaller than the IFFT size, e.g. 2048, it is possible 
to reduce number of points used in the IFFT to a number a little larger than 
the tap length. Here, such a simplified IFFT method according to the 
invention is proposed. 

Consider an IFFT filter to interpolate the channel sample while transforming 
the frequency signal into time domain. In this case, so long as the frequency 
spacing of the CTF samples is smaller than the channel coherent bandwidth, 
we do not need to know the channel information of any frequencies but the 
required number of CTF samples. Therefore, by omitting consideration of 
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such frequencies, the IFFT size and hence hardware complexity can be 
reduced. 



In an ISDB-T system, for example, the number of carriers Nc = 1405, the 
5 carrier spacing is 1/ Ts, the number of IFFT points or IFFT size, AZ/fft = 2048. 
The pilot symbols are inserted every 12 carriers, i.e. the pilot carrier spacing 
is 12/ 7s. Because the positions of the pilot symbols change periodically 
every 4 OFDM symbols, the pilot carrier spacing can be treated as 3/ Ts. 
Since the sampling interval of the incoming signal for FFT is 7c, the FIR filter 
io tap interval should be 7c too. Conventionally, after interpolation of Nc carriers 
using estimation of pilot carriers, IFFT (1405,2048) is processed. Then the 
first 128 of 2048 IFFT output samples are extracted as tap coefficients of FIR 
filter. By contrast in the present invention we hope to reduce the IFFT 
computation quantity to 256 points, i.e. extracting the tap weights from a 256 
15 IFFT output sequence. The symbol duration, 7s = 20487c, while the duration 
of 256 IFFT samples, 7; = 2567;. = 7;/8. Over frequency domain, the CTF 
samples for simplified IFFT computation with the frequency spacing of i/r 
can be obtained through interpolation of the known pilot carriers. With the 
CTF samples of virtual A/c / 8 = 176 carriers, IFFT (176,256) is processed. 
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Simulation and Results 

We now present simulations of the embodiments discussed above, and 
comparisons with simulations of the known algorithms. In the simulations, 
signals from respective sections were handled as the complex base band 
signal of equivalent low-pass systems and the calculations was conducted in 
double precision real numbers. Frequency synchronization and clock timing 
are assumed to be ideal. The relay station system is considered in a 
multipath time-invariant channel. In Table 1, the parameters used for 
simulating the coupling cancellation unit 9 in a two-path multipath channel 
are shown. The simulations shown in the figures assume that C(cd) is 
constant for all frequencies co. 



Table 1 Parameters of the SFN Repeater System Simulations 



Amplitude ratio of coupling wave : 
DU 


2dB 


Delay time of coupling wave wave : 


2.09)as 


Phase of coupling wave : 0c 


45 degrees 


Device noises of master station : 

C/Na 


40dB 


Relay station noise : C/Nb 


40dB 


Delay time of multipath : Xm 


0.984^s 



Phase of multipath : 9m 


180 degrees 


FIR filter tap length : M 


128 



Let the received signal of the receiver be ).i(/, /c), and the decision value of the 
received signal k) be d{i, /c), where k is the carrier number. Decision error 
signal err{i, k) is given as the following equation. 

err{i,k) = d{i,k)- jui^k) (30) 

The signal/noise ratio (SNR) VsnrO) of /th synnbol is defined as equation (26). 



rsvv/^(0 = 10log,o 



a: -I 



k-O 



^ err{i, k) • err * (i,k) 



k-O 



(31) 



The system performance evaluation is based on the SNR verse number of 
iteration. 

We will now explain the simplified IFFT used in the invention. Figure 6 shows 
10 the computer simulation results of the 64QAM-OFDM coupling cancellation 
unit in a multipath channel with an IFFT size of 2048 points and 256 points 
for different coupling delay times. Obviously, the cancellation unit system 
performance can be maintained even if the IFFT size is reduced from 2048 to 
256 provided that the delay time of the coupling waves is smaller than FIR 
15 filter tap length. This simulation is performed using a conventional 
cancellation unit, without CPE. 

Table 2 IFFT Computation Quantity 
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IFFT Size 

NiFFT 


Multiplication 

(N|FFT/2)l0g2N,FFT 


Addition 

N|FFTl0g2N|FFT 


2048 


11264 


22528 


256 


1024 


2048 



Table 2 shows the computation quantity of IFFT for 2048 points and 256 
points, respectively. As can be seen, the number of computations required by 
the IFFT is reduced up to 10 times after changing IFFT size from 2048 to 
256. That means the hardware complexity will be reduced greatly and 
5 processing speed will be accelerated. This method can also be applied to 
other receiver systems which use an IFFT transform e.g. decision feed back 
equalization. 

We will now discuss the application of the invention to QAM Modulation with 
scattered pilot symbols. Table 3 shows the parameters used in the simulation 
10 of a 64QAM-OFDM modem. 



Table 3 Parameters of OFDM Modem 



Modulation system 


64QAM-OFDM 


Transmission bandwidth : BW 


5.575. ..MHz 


Total number of carriers : K 


1405/Symbol 


Number of data carriers : Kd 


1287/Symbol 


Pilot carrier, etc. : Kp 


118/Symbol 


FFT size ; Nfft 


2048 


FFT sample frequency : fs 


8. 12693.. .MHz 


Effective symbol length : Ts 


252|as 


Guard intervals : Tg 


63^s(Ts /4) 


SP interpolation carrier filter 


Linear approximation 


SP interpolation symbol filter 


Same value interpolated 
within 4 symbols 


Updating rate of FIR filter 
coefficient 


4 symbols/one updating 


Error correction 


None 



Figures 7 and 8 show the results of the computer simulation of the 64QAM- 
OFDM coupling cancellation unit convergence characteristics in a multipath 
channel using the conventional method and the new proposed method 
(referred to here as Channel Pre-Estimation, CPE) under the same channel 
condition which is the working limit of the cancellation unit which uses the 



conventional nnethod. The simulation results show that the new method, 
CPE, accelerates the system convergence rate. 

Figure 9 shows the system performance using new method under its working 
limit. It is shown that the new method can improve the system performance 
greatly. The cancellation unit using CPE can work well at DUm=3dB for |_i=1 
and DUm=2dB for |_i=0.5, while conventional methods are only suitable for 
DUm=9dB and DUm=5dB respectively. 

We turn now to the application of the invention to the case of differential 
modulation. Table 4 shows the parameters used in the simulation of DQPSK- 
OFDM modem. 



Table 4 Parameters of OFDM Modem 



Modulation system 


k/4 Shift DQPSK- 
OFDM 


Transmission bandwidth 
BW 


5 575 MHz 


Total number of carriers : K 


1405/Symbol 


Number of data carriers : Kd 


1405/Syinbol 


FFT size : Nfft 


2048 


FFT sample frequency : fs 


8. 12693.. .MHz 


Effective symbol length : Tu 


252|as 


Guard intervals : Tg 


63ns(Tu/4) 


Error correction 


None 



For assessing convergent characteristics, the signal-to-noise ratio of 
respective symbols of the receiver was used after delayed detection. Figure 
10 shows the block diagram of the method of calculating the signal/noise 
ratio VsNR at the receiver. The same definition as equation (31) determines 
the signal/noise ratio Fs/vMO for /-th symbol. 

Figures 11, 12 and 13 show the result of computer simulations of 
convergence characteristics of a DQPSK-OFDM coupling cancellation unit. 
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The vertical axis represents mean SN ratio defined according to equation 
(31), and the horizontal axis is the number of iterations of the filter 
coefficients. Three cases are simulated: 

1. Using the PSK method to estimate the desired signal plus the 
conventional algorithm to cancel echoes using Equation (9). Thus, 
the methods in this case are totally conventional. 

2. Using the DPSK method to estimate the desired signal plus the 
conventional algorithm to cancel echoes using Equation (9). 

3. Using the DPSK method to estimate the desired signal plus the 
Channel Pre-estimation (CPE) method to cancel echoes using 
Equation (15). 

From Figure 11, 12 and 13, we can see that the DQPSK method is much 
better than the PSK method. This is because the conventional method 
cannot really estimate the desired signal in the multipath channel without a 
signal that consists of the channel phase. In addition, we can see that our 
most preferred method (case 3, which incorporates both of our proposed 
improvements) can achieve DUm=5dB and outperform the other methods. 
Note as DUm decreases, the system performance becomes worse. 

In summary, the presently proposed method, Channel Pre-estimation (CPE), 
improves the system performance greatly and outperforms the conventional 
method. By using channel pre-estimation, the relay station can function in an 
environment which includes a multipath channel which has deep fading. 



Moreover, the proposed simplified IFFT method reduces the required 
quantity of computations and hardware complexity. The simulation results 
show that these savings are made without reducing the performance of the 
system. 
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Claims 

1. A method of modifying the values of a plurality of digital filter 
coefficients which are used by a digital filter which is a component of a relay 
station, the relay station successively receiving a first signal, modifying the 
first signal using the digital filter to form a second signal, and transmitting the 
second signal with amplification, said first signal including a master signal 
subject to an interference and a coupling component derived from the second 
signal, the method comprising: 

(i) interrupting said transmission of the second signal; 

(ii) deriving an estimate of said interference and storing the derived 
interference estimate; 

(iii) resuming said transmission of the second signal; and 

(iv) modifying said digital filter coefficients using said stored 
interference estimate and said second signal, to reduce the difference 
between said second signal and said master signal. 

2. A method according to claim 1 in which said step (iv) is repeated 
iteratively, whereby the difference between said second signal and said 
master signal is successively reduced. 

3. A method according to claim 2 in which, upon performing each step 
(iv), a first parameter is derived characterizing the difference between said 
second signal and said master signal, and step (iv) is repeated until the value 
of said first parameter traverses a predefined convergence value. 
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4. A method according to claim 3 in which, upon said first parameter 
traversing said predefined convergence value, at least one step is performed, 
using said received signal but not using said stored interference estimate, of 
updating said digital filter coefficients to reduce the difference between said 
second signal and said master signal. 

5. A method according to claim 3 or claim 4 in which, after said first 
parameter traverses the predefined convergence value, a second parameter 
characterizing the difference between the second signal and the master 
signal is repeatedly derived, and, upon said second parameter traversing a 
predefined divergence value, said steps (i) to (iv) are repeated. 

6. A method according to any preceding claim in which the relay station 
modifies the first signal by passing it through a coupling cancellation unit 
having an input receiving the first signal and an output outputting said second 
signal, and including said digital filter and a subtracter, the input of the 
coupling cancellation unit being an input to the subtracter and an output of 
the subtracter being passed to the output of the cancellation unit and to an 
input of the digital filter, the subtracter subtracting the output of the digital 
filter from the received signal. 

7. A method according to claim 6 in which the output of the coupling 
cancellation unit includes a switch, the second signal being selectively 
transmitted via the switch to an amplifier which amplifies the second signal 
prior to transmission, said interruption of the transmission being performed by 
opening the switch. 



8. A method according to any preceding clainn in which said step of 
modifying said digital filter coefficients includes: deriving a discrete Fourier 
transform of the second signal; employing said discrete Fourier transform of 
the second signal and a discrete Fourier transform of said estimated 
interference to derive correction values for said digital filter coefficients 
expressed in the frequency domain; and transforming said correction values 
into the time domain by an inverse Fourier transform operation, the digital 
filter coefficients being modified based on the time domain correction values. 

9. A method according to claim 8 in which each of the digital coefficients 
of the digital filter modifies the first signal at pairwise time delay intervals Tp, 
the product of Tp and the number of digital coefficients being substantially 
equal to a maximum delay time of a mulitpath interference component of the 
first signal. 

10. A method according claim 9 in which the inverse Fourier transform 
generates an number of output time domain values equal to or slightly larger 
than the number of digital coefficients. 

11. A method according to claim 9 or claim 10 in which the inverse Fourier 
transform operates on a number of frequency domain values lower than the 
number of subcarriers used in practice. 

12. A method according to any preceding claim in which said master 
signal is in a format in which symbols are coded by differences between 
consecutive segments of the master signal, and said adjustment of the digital 
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filter coefficients minimizes an estimated error in the estimation of differences 
between consecutive segments of the master signal. 

13. A method according to claim 12 in which the format of the master 
signal encodes symbols by phase differences between consecutive 

5 segments. 

14. A method of modifying the values of a plurality of digital filter 
coefficients for use by a digital filter which is a component of a relay station, 
the relay station successively receiving a first signal, modifying the first signal 
using the digital filter to form a second signal, and transmitting the second 

10 signal with amplification, said first signal including a master signal subject to 
an interference having a maximum delay time, each of the digital coefficients 
of the digital filter modifying the first signal at pairwise time delay intervals Tp, 
the method including: 

deriving a discrete Fourier transform of the second signal; 

15 employing said discrete Fourier transform of the second signal to 

derive correction values for said digital filter coefficients expressed in the 
frequency domain; and 

transforming said correction values into the time domain by an inverse 
Fourier transform operation, the digital filter coefficients being modified based 
20 on the time domain correction values, 

wherein: 
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the product of Tpand the number of digital coefficients is substantially 
equal to the maximum delay time of the interference component of the first 
signal, and 

the inverse Fourier transform generates an number of output time 
domain values equal to or slightly larger than the number of digital 
coefficients. 

15. A method of modifying the values of a plurality of digital filter 
coefficients for use by a digital filter which is a component of a relay station, 
the relay station successively receiving a first signal, modifying the first signal 
using the digital filter to form a second signal, and transmitting the second 
signal with amplification, said first signal including a master signal and a 
coupling component, the master signal being in a format in which symbols 
are coded by differences between consecutive segments of the master 
signal, and the coupling component being derived from the second signal, 
the method comprising: 

(i) deriving an estimate of differences between consecutive 
segments of the master signal; 

(ii) estimating an error in said estimation of differences between 
consecutive segments of the master signal; and 

(iii) modifying said digital filter coefficients to reduce said estimated 
error. 



16. A coupling cancellation device for a relay station, the relay station 
receiving a first signal and transmitting it to a coupling cancellation device, 
receiving a second signal fronn the coupling cancellation device, transmitting 
5 the second signal to an amplifier and transmitting the amplified second 
signal, 

the coupling cancellation device including a digital filter based on a 
plurality of filter coefficients and arranged to generate a correction signal for 
combination with the first signal for reducing interference in the second 
10 signal, a digital signal processing unit which determines the filter coefficients, 
and a switch which selectively transmits the second signal to the amplifier of 
the relay station; 

whereby the switch can be controlled not to transmit the second signal 
when it is desired to estimate multipath interference in the first signal. 

15 17. A coupling cancellation device according to claim 16 in which the 
digital filter receives and the digital signal processing unit each receive said 
second signal, the coupling cancellation device further including a subtraction 
circuit receiving the correction signal from the digital filter and subtracting the 
correction signal from the first signal to generate the second signal. 

20 18. A coupling cancellation device according to claim 17 in which the 
digital signal processing unit transforms the second signal into the frequency 
domain using a FFT to extract signals over pilot carriers for estimation of the 
channel transfer function, the digital signal processing unit including a switch 
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for selectively transmitting the channel transfer function to a multipath 
interference estimation unit or to a digital filter coefficient determination unit. 
19. A coupling cancellation device for a relay station, the relay station 
receiving a first signal and transmitting it to a coupling cancellation device, 
5 receiving a second signal from the coupling cancellation device, transmitting 
the second signal to an amplifier and transmitting the amplified second 
signal, the first signal being subject to an interference having a maximum 
delay time, 

the coupling cancellation device including: 
10 a digital filter based on a plurality of filter coefficients and arranged to 

generate a correction signal for combination with the first signal for reducing 
interference in the second signal, each of the digital coefficients of the digital 
filter modifying the first signal at pairwise time delay intervals Tp, and 

a digital signal processing unit which determines the filter coefficients 

15 by: 

deriving a discrete Fourier transform of the second signal; 

employing said discrete Fourier transform of the second signal to 
derive correction values for said digital filter coefficients expressed in the 
frequency domain; and 

20 transforming said correction values into the time domain by an inverse 

Fourier transform operation, the digital filter coefficients being modified based 
on the time domain correction values, 

wherein: 
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the product of Tp and the number of digital coefficients is substantially 
equal to the maxinnum delay time of the interference component of the first 
signal, and 

the inverse Fourier transform generates an number of output time 
5 domain values equal to or slightly larger than the number of digital 
coefficients. 

20. A coupling cancellation device for a relay station, the relay station 
receiving a first signal and transmitting it to a coupling cancellation device, 
receiving a second signal from the coupling cancellation device, transmitting 
10 the second signal to an amplifier and transmitting the amplified second 
signal, the master signal being in a format in which symbols are coded by 
differences between consecutive segments of the master signal, 

the coupling cancellation device including a digital filter based on a 
plurality of filter coefficients and arranged to generate a correction signal for 
15 combination with the first signal for reducing interference in the second 
signal, and a digital signal processing unit which determines the filter 
coefficients, 

the digital signal processing unit being arranged to derive an estimate 
of differences between consecutive segments of the master signal, estimate 
20 an error in said estimation of differences between consecutive segments of 
the master signal, and modify said digital filter coefficients to reduce said 
estimated error. 
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21. A relay station incorporating a coupling cancellation unit 
according to any of claims 16 to 20. 
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Abstract 

Methods and Apparatus for Reducing Signal Degradation 

A single frequency relay station receives a signal from a master station, and 
retransmits it. Before retransmitting the signal, the relay station uses a digital 
5 filter to subtract from the signal components which arise from multipath 
interference and coupling interference. The coefficients used by the digital 
filter are derived from characteristics of the multipath interference and the 
overall transfer function of the relay station. These are derived by turning off 
the retransmission, so that the multipath interference can be estimated from 
10 the received signal, and commencing the retransmission again, to determine 
the transfer function. A simplified inverse fast Fourier transform is used to 
simplify the calculations. 
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